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ABSTRACT 

The  angular  distributions  of  78  +  5  Mev  i t~  and  ir 

mesons  scattered  elastically  from  copper  have  been  measured. 

The  experimental  results  are  compared  with  the  distributions 

found  from  a  range  of  exact  optical  model  phase  shift 

calculations  using  a  complex  attractive  square  well  potential 

1/3  —13 

for  r  <(  Rq  (where  Rq  =  1.4A  x  10  cm)  and  the  Coulomb 
potential  for  r  )>  Rq.  The  effect  of  adding  to  the  potential 
a  term  depending  upon  the  gradient  of  the  nuclear  density, 
and  thereby  introducing  a  discontinuity  in  the  matching  of 
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logarithmic  derivatives  at  r  =  RQ,  is  also  explored.  The 
optical  model  as  well  as  several  modified  Born  approximation 
calculations  show  reasonable  agreement  with  the  experimental 
curve  in  the  forward  direction  and  in  the  region  of  the  first 
diffraction  maximum.  At  backward  angles  the  calculated  curves 
show  pronounced  maxima  and  minima  which  do  not  appear  experi¬ 
mentally,  perhaps  because  incoherent  transitions  at  these 
angles  mask  the  true  coherent  distribution  in  the  experimental 


measurement. 
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I.  INTRODUCTION 

1  + 

The  elastic  scattering  of  78  +  5  Mev  i r  and  tr  by  copper 


The  sum  of  true  elastic  plus  nearly  elastic  scattering  is 
actually  measured. 


was  measured  to  examine  the  pion  angular  distribution  produced 

by  a  medium  -  heavy  nucleus.  This  is  one  part  of  a  program  to 

explore  variations  in  positive  and  negative  pion  differential 

cross  sections  over  the  range  of  available  target  nuclei.  A 

2 

recent  measurement  of  elastic  pion  scattering  from  lithium 
revealed  a  pronounced  minimum  and  a  strong  backward  rise 


R.  E.  Williams,  J.  Rainwater,  A.  Pevsner,  Phys.  Rev.  101, 
412,  (1956) 


characteristic  of  the  elementary  pion-single  nucleon  angular 
distributions.  For  the  (Z  =  3)  lithium  nucleus  an  analysis 
was  made  in  terms  of  a  coherent  addition  of  elementary 
scattering  amplitudes  using  a  modified  impulse  approximation 

3 

analysis.  An  earlier  experiment  had  shown  no  pronounced 


A.  Pevsner,  J.  Rainwater,  R.  E.  Williams,  and  S.  J. 
Lindenbaum,  Phys.  Rev.  100,  1419  (1955) 


minimum  and  only  a  small  backward  rise  for  elastic  pion 
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scattering  from  the  somewhat  heavier  (Z  =  13)  aluminum  nucleus. 

An  optical  model  phase  shift  calculation  using  a  complex  square 
well  potential  showed  strong  maxima  and  minima  for  aluminum  which 
did  not  appear  in  the  experimental  curves,  but  otherwise  the 
general  characteristics  of  the  calculated  curves  at  angles  below 
that  of  the  first  diffraction  minimum  could  be  approximately 
matched  to  the  experimental  results.  For  the  present  experiment, 
copper  (Z  =  29)  was  chosen  as  a  material  of  medium  Z  whose  nuclear 
radius  approximately  equals  the  pion  mean  free  path  for  absorption 
in  nuclear  matter.  For  the  heavier  copper  nucleus  the  complex 
square  well  might  be  expected  to  better  represent  the  overlapping 
fields  of  the  individual  nucleons,  and  the  experiment  was  carried 
out  with  relatively  high  angular  resolution  and  good  statistical 
accuracy  in  an  attempt  to  avoid  masking  whatever  optical  model 
diffraction  effects  might  be  present.  The  measured  angular 
distribution  shows  some  suggestion  of  a  diffraction-like  oscil¬ 
lation  effect  in  the  region  of  40°  to  80°.  Optical  model  calcu¬ 
lations  have  been  made  using  a  complex  square  well  potential. 

The  effect  of  including  an  added  term  which  depends  upon  the 

4 

gradient  of  the  nuclear  density  was  also  tested.  When  a  square 


L.  S.  Kisslinger,  Phys.  Rev.  98,  761  (1955) 


well  potential  is  used  the  gradient  term  alters  the  boundary 
condition  for  matching  inside  and  outside  solutions  at  the  nuclear 
surface . 
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II.  EXPERIMENTAL  ARRANGEMENT 

Thp-  scintillation  counters  and  coincidence  circuits  used 

to  measure  pion  angular  distributions  from  the  copper  target 

are  the  same  ones  previously  described  for  the  pion  experiments 

with  lithium  and  aluminum  targets.  The  arrangement  of  the 

scintillation  counter  telescopes  and  a  description  of  the 

2  3 

electronic  circuitry  have  been  given  in  earlier  papers  ’  . 

In  brief  summary,  the  telescopes  are  place!  behind  a  focusing 

magnet  in  the  80  Mev  pion  exit  channel  of  the  380  Mev  Nevis 

synchrocyclotron,  as  shown  in  Figure  1.  To  measure  angular 

distributions,  counters  3  and  4  can  be  rotated  vertically  about 

—8 

the  target  position  as  a  pivot  point.  Three  fast  (10  sec) 

coincidence  circuits  are  used  between  various  counter  pair 

combinations  in  the  two  telescopes  in  such  a  manner  as  to 

reduce  background  due  to  random  coincidences  between  counters 

~7 

to  a  low  level;  while  a  slower  (10  sec)  triple  coincidence, 
operating  on  the  output  of  the  faster  circuits,  establishes 
the  overall  coincidence  rate. 

III.  EXPERIMENTAL  RESULTS 

The  measured  angular  distributions,  together  with  the 
spreads  in  angular  resolution,  are  shown  for  78  Mev  negative 
pions  in  Table  I,  and  for  78  Mev  positive  pions  in  Table  II. 

The  +  quantities  listed  with  the  cross  sections  are  the 
statistical  standard  deviations.  The  results  are  also  shown 
in  Figure  2.  The  many  corrections  to  the  raw  counting  data 
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required  to  arrive  at  the  differential  cross  sections  have  been 

2  3 

J  J  4  «  r-,  4-  n  4  ~l  A  n  no  r»1  4  qy»  r\o  nor»e  *  TjlGSG  iHClud.6 


described  in  detail  in  the  earlier  papers 


adjustment  for  counting  efficiency  in  the  electronics,  multiple 
Coulomb  scattering  in  the  target  and  copper  absorber,  incident 
beam  width,  finite  counter  size,  ir-p.  decay  after  traversing  the 
target,  electron  and  p  meson  contamination  in  the  incident  beam, 
stopping  power  of  the  target,  nuclear  absorption  of  full  energy 
mesons  in  the  target  and  copper  absorber,  low  energy  pions  in 
the  incident  beam,  etc.  Plots  of  the  values  included  in  Tables  I 
and  II  are  shown  in  Figure  2.  The  ir~  curve  shows  some  evidence 
of  a  diffraction-like  oscillation  effect  near  the  bottom  of  the 


initial  decrease  of 


in  going  from  small  to  moderate  angles. 


The  points  near  60  appear  as  a  relative  bump  on  the  curve  that 
would  be  constructed  by  interpolating  from  the  points  on  either 
side.  This  is  the  first  evidence  of  such  diffraction  * ’structure 
in  the  scattering  of  pions  from  complex  nuclei,  and  provides  an 
interesting  feature  for  comparison  with  the  results  of  optical 
model  phase  shift  calculations. 


!V.  NUCLEAR  MODELS 

5  6 

Watson  and  Francis  have  given  general  arguments  for  the 


5  K.  M.  Watson,  Phys.  Rev.  8-9,  575  (1953) 

6  N.  C.  Francis  and  K.  M.  Watson,  Phys.  Rev.  92,  291  (1953) 


validity  of  nuclear  optical  model  potentials  for  pion  scattering 
from  nuclei.  It  is  expected  for  light  nuclei  that  the  results 
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will  be  strongly  influenced  by  characteristics  of  the  elementary 

4 

pion-nucleon  amplitudes  ,  whereas  for  heavier  nuclei  an  optical 

model  analysis  using  a  complex  interaction  potential  should 

2 

become  more  valid.  In  a  preceding  paper  of  this  series  a 
specific  nuclear  model  was  used  with  the  impulse  approximation 
to  analyze  the  results  of  pion  scattering  from  lithium  in  terms 
of  the  angularly  dependent  elementary  pion-nucleon  scattering 
amplitudes.  To  obtain  definite  individual  particle  nucleon 
wave  functions,  a  simple  harmonic  nuclear  well  was  assumed  for 
the  nucleons.  The  individual  terms  in  the  coherent  addition 
of  elementary  scattering  amplitudes  from  the  various  nucleons 
were  weighted  according  to  the  ability  of  the  struck  nucleon 
to  absorb  the  momentum  recoil  and  remain  in  the  same  state  in 
the  simple  harmonic  nuclear  well.  Several  obvious  adjustments, 
such  as  relativistic  solid  angle  and  phase  space  corrections, 
were  made.  Corrections  were  also  made  at  backward  angles  for 
the  momentum  distribution  of  the  nucleons.  By  also  taking 
into  account  the  energy  dependence  of  the  elementary  f(0), 
fairly  good  agreement  with  the  experimental  angular  distribution 
was  obtained.  The  theoretical  calculation  was  quite  sensitive 
to  the  choice  of  such  parameters  as  nuclear  size  and  the 
approximations  used  in  applying  the  correction  factors.  A 
best  match  with  the  experimental  data  corresponded  to  ’ ’reasonable * 
choices  of  these  parameters. 

Kisslinger  has  constructed  nuclear  scattering  potentials 

4 

directly  from  the  elementary  s  and  p  wave  amplitudes  and  has 
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shown  the  importance  of  a  surface  term  which  arises  from  the  p 

wave  contribution  to  the  elementary  f(9).  The  magnitude  of  the 

surface  term  depends  upon  S'  p,  where  p(r)  is  the  nuclear  density 

distribution.  Variations  in  both  range  and  shape  of  the  nuclear 

7 

density  distribution  about  various  favored  distributions  have 


7  D.  L.  Hill  and  K.  W.  Ford, 

Phys. 

Rev.  94,  1617  (1954) 

8 

been  made  to  fit  both  the  experimental  p  mesic  x-ray  data  and  the 

g 

V.  Fitch  and  J.  Rainwater, 

Phys. 

Rev.  92,  789  (1953) 

9  10 

electron  scattering  data.  } 

Such 

studies  have  led  to  a  favored 

9 

Yennie ,  Ravenhall,  and  Wilson,  Phys.  Rev.  95,  500  (1954) 

10  R.  Hofstadter,  Proceedings  of  the  Fifth  Annual  Rochester 
Conference  on  High-Energy  Physics,  Interscience  Publishers, 
Inc.,  New  York,  1955. 


nuclear  density  distribution  which  is  relatively  uniform  in  the 
center  region  of  the  nucleus  but  tapers  gradually  to  zero  at  the 
nuclear  ftsurface'’.  It  should  be  noted,  however,  that  the  fx  mesic 
x-ray  and  the  electron  scattering  experiments  give  the  charge 
distribution  in  the  nucleus,  whereas  here  we  are  interested  in  the 
average  pion  interaction  as  a  function  of  r.  This  may  be  non¬ 
linear  in  nucleon  density,  may  have  a  range  extension,  and  may 
reflect  differences  in  the  neutron  and  proton  distributions. 
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V.  THEORY  AND  CALCULATIONS 

We  have  performed  phase  shift  optical  model  calculations 
using  a  constant  complex  square  well  potential  within  the 
nuclear  radius  (r  -  RQ  =  1 .4A1// x  lO-13  cm.),  a  Kisslinger 
type  term  at  the  surface,  and  a  Coulomb  potential  beyond  the 
nuclear  radius.  Exact  solutions  are  obtained  to  the  Klein 
Gordon  equation  in  the  form 


v2f  + 


(e  -  vy 


2  4 
m  c 


•^  =  0 


Inside  the  nuclear  surface  V  =  Y^  +  1Y2,  at  the  surface 

V  -ft  =  -BV  p  •  L?  ip  where  B  is  a  negative  constant,  and  beyond 

-Ze2 

the  surface  V  =  — —  .  A  range  of  inside  potentials  was  chosen 
taking  into  account  the  increased  Coulomb  and  nuclear  contri¬ 
butions  from  the  heavier  copper  nuclei  compared  with  aluminum, 

for  which  the  v~~  potentials  V-,  =  -30  Mev  and  V,-  -  -10  to  -25 

1 

3 

Mev  seemed  roughly  appropriate.  The  following  sets  of  values 
were  used: 


(a) 

V1  ~ 

-25 

Mev, 

V2  = 

-20 

Mev 

(b) 

V1  = 

-35 

Mev, 

V2  " 

-20 

Mev 

<c) 

V1  " 

-45 

Mev, 

V2  " 

-20 

Mev 

(d) 

V1  " 

-25 

Mev, 

V2  = 

-15 

Mev 

Nuclear  phase  shifts  were  found  in  the  conventional  way 
by  matching  logarithmic  derivatives  of  the  regular  and  irregular 
Coulomb  solutions  (for  r  )>  Rq)  to  the  spherical  Bessel  functions 
of  complex  argument  inside  the  nucleus.  The  surface  potential 
affected  only  the  matching  of  logarithmic  derivatives  of  the 
wave  functions  at  the  nuclear  boundary.  To  be  sure  that 
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significant  partial  wave  contributions  were  not  neglected,  ft 
values  from  0  to  6  were  included  in  the  calculations.  Actually 
contributions  from  partial  waves  with  ft  )>  5  should  be  small,  and 
the  results  showed  negligible  contribution  for  ft  )>  6. 

For  the  radial  part  of  the  inside  solution  the  complex 
propagation  constant  is  given  by 


2 

1  "  eV 


2  4 
m  c 


(2) 


where  V  =  V,  +  iVQ  has  been  assumed  somewhat  arbitrarily  to 
transform  as  the  fourth  component  of  a  four  vector,  in  analogy 
to  the  Coulomb  potential  for  the  outside  solution.  Wave  functions 
for  the  inside  Schrodinger  equation  are  then  the  spherical  Bessel 
functions 


h 

r 


J  (hr) 
i+£  1 


(3) 


from  which,  by  the  usual  recurrence  relations,  the  logarithmic 
derivatives  at  r  =  R^  were  found  for  ft  values  from  0  to  6. 

If  the  small  quadratic  V  term  is  neglected  in  (1),  the 
outside  Schrodinger  equation  reduces  to  the  usual  form 

^7  ifj  +  (k0  -  _ 0  )  ip  =  0 


where 


and 


k„2  -  2s  E<  (1  +  -14  ) 


0  "? 


2mc 


E '  =  E  -  me  =  78  Mev 
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T  /*?  •? 

Choosing  Rq  =  1.414A  x  10  cm.,  then  =  4.8  and 


a  = 


mZ 


2 


h^k 


(1  + 


E' 


me 


)  =  -.27487 


for  78  Mev  it"  mesons  incident  upon  a  copper  nucleus.  In  order 
to  represent  an  outgoing  scattered  wave  superimposed  upon  an 
incident  wave  unmodified  by  the  nuclear  potential,  the  outside 
wave  function  must  have  the  form 


yp  i<5  p 

c£*^U,k0R0)  =  F^(a,k0R0)  +  [G^ (a , kQR0>  +  iF^ (a, kQR0)  ]  e  sinb 

(5) 

where  F^  and  G£  are  the  regular  and  irregular  hypergeometric 
wave  functions.  Values  of  F  £  were  indirectly  obtained  from 
the  NBS  tables11  which  list  <3?^(a,k0R0>  and  CQ(a),  related  to 

11  Tables  of  Coulomb  Wave  Functions,  Volume  I,  National 

Bureau  of  Standards,  Applied  Mathematics  Series  17,  1952. 


by 


and 


F£  =  C£(a)  (k0R0) 


£+1 


C £(a)  = 


L  2  2 

y£  +  a 

£(2£  +  1) 


C£(-a)  =  em  C£  (a) 


(a,k0R0) 

Cf-1 
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12 

Values  of  are  not  tabulated,  except  for  positive  values  of  a. 

12  I.  Bloch,  M.  M.  Hull,  Jr.  A.  A.  Broyles,  W.  G.  Bouricius, 

B.  E.  Freeman,  and  G.  Breit,  Rev.  Mod.  Phys.  23,  147  (1951) 


Since  our  a  =  -.27487,  the  existing  tables  could  not  be  extrapolated 
from  a  =  0  with  sufficient  accuracy.  Therefore  G^a, kQR0)  was 
calculated  from  the  very  lengthy  and  slowly  converging  series 
expansions  for  G^.  Other  G^  values  were  then  found  from 


„  HI2  *  a2)-i  .  Vi- 1 

G  =  -  +  - 

X*  m  in 


l-l 


£-1 


F*,  was  calculated  with  the  use  of  the  recurrence  relations 


£F^  =  Vri2+a2  -  ( -  +  a)Fj>  for  £  ^  0 


k0R0 


(i+l)F*  =  [ 


_  r  (i+lV 


koRo 


+  a]  F£  -  [(f+1)2  +  c?Y  V&+1  for  £  -  0 


G |  can  then  be  found  from  the  Wronskian  equation 


riGi  -  FiG]  "  1 


13 

The  equation  for  the  scattering  amplitude  is 


L.  I.  Schiff,  Quantum  Mechanics,  sec.  20,  (McGraw  Hill,  New 
York,  1949) 
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oo  1  i  (2ct  n+6  /)) 

f(0)  =  f  (0)  +  2  i  (24+1)  e  *  sin6  »  P.  (cos0)  (6) 

c  i«0  *  £ 


where 


fc(e> 


a 


2  g 

-ia  In  (sin  ^ )  +  2iaQ 


2ksin 


2” 0 


11 

Values  for  can  be  had  from  the  tabulated  values  by 

means  of  the  recurrence  relation 

—1  a 

°Ji  =  a£~l  +  tan  ^ 

If  the  nucleon  density  is  assumed  constant  within  the 

nuclear  radius  (corresponding  to  a  square  well  model)  but 

drops  to  zero  within  a  region  y  -*■  0  about  the  assumed  radius, 

4 

the  matching  condition  at  the  boundary  becomes 


dr 


(7) 


2E 


where  n  =  1  +  2  Bp  =  0.473 


K  c 


meson  total  energy 

nucleon  density  normalized  so  / pdV  =  1 


E 
P 

B  =  (2i r)3  At 
A 
t 


total  number  of  nucleons 
3 


2 

-E  c 


(2ir)  E  k3 


\+  'C  E3/S!  ^  ^  T3/2  ^ 


J3/2 


3/2 


-.3611  x  10  ^erg  cm3 
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E3/2 

T3/2 


e  ++sin5  for  both  isotopic  spin  =  3/2  and 

4~r 

j  =  3/2  in  pion-nucleon  scattering 
=  spin  projection  operator  for  j  =  3/2  in  the  pion- 
nucleon  system 

=  isotopic  spin  projection  operator  for  the  T  =  3/2 
state  of  the  pion-nucleon  system. 


Thus  the  new  boundary  matching  condition,  taking  account  of  the 
logarithmic  discontinuity,  is 


=  (!  -  n)  + 

l  k0R0 


n 


k0 


(8) 


In  our  calculation  we  have  used  n  =  0.5,  0.75,  and  1.00  instead 
of  the  calculated  0.473  for  the  following  reasons: 

(a)  Several  choices  of  n  values  will  better  demonstrate 
the  effect  of  the  surface  term. 

(b)  An  exact  choice  for  n  is  unnecessary,  since  (7) 
above  has  not  been  corrected  for  the  difference  in  available 
phase  space  between  pion-nucleon  and  pion-nucleus  center  of 
mass  systems,  or  for  the  off-energy  shell  scattering  which  is 
actually  taking  place.  Furthermore,  the  p  state  interaction 
between  pion  and  nucleon  gives  rise  to  a  mass  like  term  and 
an  ordinary  potential  term  in  the  pion-nucleus  Hamiltonian, 
in  addition  to  the  vp  term.  A  consistent  treatment  would 
include  expressions  for  these  two  additional  terms;  we  have 
done  so  only  through  our  phenomenological  determination  of 
the  complex  square  well  depths. 
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Nevertheless,  a  change  in  the  angular  distribution 
resulting  from  a  change  in  n  will  at  least  be  indicative  of 
the  effect  of  the  surface  potential  contribution  to  the 
scattering. 

The  equation  determining  the  nuclear  phase  shifts  6^ 
from  matching  at  the  boundary  is 

e2i6<  ,  +  1Gi>  -  +  iGi>  hi  (9) 

-(F)  -  iG))  +  (F^  -  iG,) 

From  (9)  5^,  and  hence  the  partial  wave  coherent  scattering 
amplitudes  (6),  can  be  found. 

VI*  results  OF  CALCULATIONS 

The  calculations,  carried  out  on  a  desk  calculator,  are 
most  easily  studied  by  plotting  f(0)  in  the  complex  plane. 
Figures  3a  and  3b  show  the  complex  f(0)  plots  for  =  -35 
Mev,  V2  =  -20  Mev,  and  boundary  matching  factors  n  =  0.50, 
0.75,  and  1.00.  The  corresponding  plot  of  the  differential 
cross  section  is  shown,  along  with  the  experimental  curve, 
in  Figure  4.  Similar  complex  plots  for  f (0)  derived  from 
V1  =  -45  Mev,  V2  =  -20  Mev  and  the  same  range  of  n  values 
are  shown  in  Figures  5a  and  5b.  The  corresponding  angular 
distributions  are  shown  in  Figure  6.  Figure  7  shows  the 
calculated  angular  distributions  for  all  four  sets  of  complex 
potential  wells  previously  mentioned,  but  for  no  discontinuity 
in  logarithmic  derivatives  at  the  nuclear  boundary  (n=l) . 
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It  is  obvious  from  Figure  7  that  the  calculated  curves 
heavily  overemphasize  the  diffraction  maxima  and  minima  in  the 
angular  distribution,  relative  to  our  experimental  curve.  However, 
it  is  difficult  to  say  how  strongly  the  experimental  curve 

excludes  such  sharp  prominences  for  strictly  coherent  scattering. 

< 

In  our  experimental  measurements,  states  of  nuclear  excitation  9 
Mev  can  contribute  to  the  measured  cross  section.  Thus  incoherent 
transitions  such  as  spin  flip  and  slightly  inelastic  scattering 
may  be  significant  factors  in  the  10  mb.  differential  cross  sections 
observed  at  backward  angles.  Furthermore,  our  angular  resolution 
of  +  3.9°  to  6.4°  would  tend  to  smooth  out  any  sharp  maxima  and 
minima  which  might  exist  in  the  purely  coherent  angular  distri¬ 
bution.  Finally,  a  considerable  damping  of  the  calculated 
diffraction  maxima  and  minima  can  probably  be  expected  from  a 
more  reasonable  nuclear  model  which  abandons  the  sharp  edged 
square  well  for  a  nuclear  density  distribution  tapered  toward 
the  edge.  Calculations  with  such  a  model  require  the  use  of 
electronic  computers  and  were  not  attempted  here. 

The  theoretical  and  experimental  curves  are  best  compared 
on  the  basis  of  the  diffraction-like  oscillation  in  the  region 
from  40°  to  80°  and  the  height  and  slope  of  the  curves  at  angles 
smaller  than  40°.  From  Figure  7  it  would  seem  that  best  agree¬ 
ment  with  the  angular  location  of  the  first  ' 'hump'*  in  the 
experimental  curve  might  be  obtained  with  the  complex  potential 

V  =  (-45  -  i  20)  Mev,  although  the  first  maximum  in  the 

V  =  (-35  -  i  20)  Mev  calculated  curve  is  only  slightly 
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•displaced  toward  larger  nnglss  rrsZLaxxvs  xd  Pne  sxpsrixnsnxa.1 
•curv-s;,  (If  21  smaller  jnidaar  raiizs  nr  values  -of  n  less  than 
unity  are  assumed,  -the  ls.ga.ti.xm  of  Phe  calculated  first 
■diffraction  maximum  mil  be  shf-fPed  toward  still  larger  -angles.. ) 
The  maximum  is  fixe  calculatsd  mrraas  at  53d  to  60°  is  net 
simply  tie  result  of  interference  befwesn  tie  Coulomb  seat  ter  lag 
amplitude  and  a  moaotogi-cally  degreasing  nuclear  scattering 
amplitude,  The  maximum  -would  still  he  present  as  part  of  tie 
diffraction  pattern  produced  by  tie  complex  spusre  well  if  the 
Coulomb  potential  were  absent,  although  it  would  then  be  some-- 
what  mociif '  *d  in  amplitude  and  shifted  in  angular  location., 

At  angles  smaller  than  40°,  both  the  ¥  =  (-45  —  i  20)  ltev 
and  ¥  •=  (-35  -  i  20)  Ker  calculated  cur-yes  agree  fairly  well 
with  the  experimental  values.  In  this  region  the  partial  wave 
scattering  amplitudes  are  adding  roughly  in  phase  and  inter¬ 
ference  effects  between  the  various  partial  wave  contributions, 

which  may  he  sensitive  to  the  assumed  nuclear  shape,  are  less 

cl  <5 

important,  faking  into  account  the  slight  lowering  of. 
brought  about  by  adding  a  surface  term  to  the  potential 
(decreasing  n  from  unity),  slightly  better  agreement  fox-  a 
range  of  n  values  at  forward  angles  is  had  with  the  potential 
V  =  (-45  -  i  20)  Mev  than  with  V  =  (-35  -  i  20)  Mev.  In  the 
backward  direction  the  curve  for  V  =  (-45  -  i  20)  Mev  also 
stands  higher  and  closer  to  the  experimental  curve,  hut  it  is 
in  this  region  that  incoherent  scattering  may  be  dominant 
and  comparisons  therefore  less  valid. 
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In  comparing  the  calculated  and  experimental  angular  distri¬ 
butions  it  is  difficult  to  assign  a  particular  choice  to  n,  which 
determines  the  discontinuity  in  logarithmic  derivatives  at  the 
nuclear  boundary.  The  choice  of  n  is  sensitive  to  the  locations  of 
maxima  and  minima  in  the  backward  direction,  but  since  these  are 
not  clearly  apparent  in  the  experimental  curve  they  are  not  useful 
criteria.  As  n  is  decreased  from  unity  the  first  minimum  and 
subsequent  maximum  in  the  calculated  curves  are  damped  and  displaced 
slightly  toward  larger  angles.  The  plots  of  f(0)  in  the  complex 
(q  +  ip)  plane  in  Figures  3  and  5  show  how  this  is  brought  about. 

The  complex  f(0)  plots  also  show  that  changes  in  location  and 
magnitude  of  maxima  and  minima  at  larger  angles  are  not  easily 
predicted  as  n  varies.  The  locations  and  depths  of  the  minima 
are  determined  by  where  and  how  close  the  f(0)  plots  swing  toward 
the  origin  as  n  is  varied,  and  from  Figures  3  and  5  it  is  seen 
that  these  features  vary  in  a  complex  manner.  In  Figures  4  and  6 
a  higher  backward  cross  section  is  obtained  for  the  smallest  n 
(n  =  0.5),  but  this  value  of  n  does  not  allow  a  good  match  with 
the  experimental  curve  in  the  region  of  the  first  diffraction 
maximum.  Values  of  n  from  0.75  to  1.00  appear  to  be  favored  when 
matching  at  the  first  diffraction  maximum. 

4* 

Optical  model  calculations  were  not  made  to  fit  the  ir 
experimental  distribution.  However,  the  statistically  weak 
experimental  points  at  50°  and  60°  suggest  that  a  maximum 
may  not  be  present  in  the  7T+  curve  in  this  region. 
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The  total  reaction  cross  section  has  also  been  calculated 

from  the  partial  wave  nuclear  phase  shifts  for  each  of  the 

real  and  imaginary  square  well  potentials  and  boundary 

matching  factors  already  discussed.  Results  are  shown  in 

Table  III.  For  most  values  of  the  potentials,  is  slightly 

larger  than  the  copper  geometric  cross  section  of  1006  mb 

1/3  -13 

(assuming  Rq  =  1.414  A  x  10~  cm).  A  similar  result  was 
found  in  the  optical  model  study  of  pion  scattering  from 

3 

aluminum  .  Although  the  geometric  cross  section  is  the  limit 

14 

for  absorption  in  the  simple  Fernbach,  Serber,  Taylor  treatment, 


14 


Fernbach,  Serber,  and  Taylor,  Phys.  Rev.  75, 


1352  (1949) 


much  larger  than  geometric  cross  sections  become  possible  with 
a  partial  wave  analysis.  The  corresponding  mean  free  paths 


for  absorption,  A.  =  ,  are  all  in  the  neighborhood  of 

a  ZKr, 

-13  2 
4  x  10  cm. 


We  have  also  made  two  modified  Born  approximation 
calculations  for  scattering  from  a  copper  target  (linear 
superposition  of  scattering  amplitude  contributions  from  the 
different  nucleons)  in  order  to  compare  with  the  results  of  the 
exact  phase  shifts  of  the  optical  model  as  well  as  to  aid  in 
the  choice  of  well  depths  for  the  lengthier  optical  model 
calculations.  Although  the  normal  first  Born  approximation 
does  not  take  into  account  attenuation  effects  of  the  difference 
in  propagation  constants  inside  and  outside  the  nucleus,  it 
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does  make  allowance  for  non-isotropic  scattering  from  each  element 
of  nuclear  matter,  and  coherently  combines  the  scattering 
contributions  from  all  parts  of  the  nucleus  according  to  a 
weighting  function  f^CS)  which  is  found  from  the  nuclear  density 
distribution.  Thus  the  total  scattering  amplitude  f (0)  can  be 
written  as  f  (0)  f, (0),  where  f  (0)  is  the  nuclear  volume  times 

3.  D  3 

the  Born  contribution  to  the  scattering  amplitude  per  unit  volume 

of  nuclear  matter,  and  f^(0),  the  nuclear  density  weighting 

function,  gives  the  effect  of  the  finite  nuclear  size.  It  is 
3  . 

given  by 

MS)  -  A2  p(r)  dr 

D  0  qr 

Our  modifications  to  this  usual  Born  approximation  have  already 
15 

been  described,  and  are  briefly  noted  below. 


15 


A.  Pevsner  and  J.  Rainwater,  Phys.  Rev.  100,  1431 


(1955) 


Modified  Born  Approximation  A 

Q 

f^(0)  is  modified  to  the  extent  of  using  =  2k^sin,| 

0 

instead  of  q„  =  2k„sin^-,  where  k^,  and  k .  =  k..  +  i  k_  are  the 

values  of  k  at  r  =  oo  and  r  <(  Rq  respectively.  This  change  brings 

the  calculated  diffraction  minima  to  about  the  same  angular 

positions  as  in  the  exact  phase  shift  calculations.  Also  the 

interior  parts  of  the  nucleus  are  weighted  less  heavily  by  using 

-k2(R0~r) 

an  attenuation  factor  e  .  Thus 
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Plots  for  the  modified  Born  approximation  A  are  shown  in 
Figure  8  for  three  of  the  four  sets  of  potentials  for  which 
the  exact  phase  shift  calculations  were  performed.  Comparison 
with  Figure  7  shows  that  the  phase  shift  and  modified  Born 
approximation  results  are  very  similar.  The  minima  are  more 
pronounced  and  the  entire  curve  somewhat  lower  in  the  Born 
case.  Angles  for  the  minima  are  nearly  the  same,  but  slightly 
smaller  in  the  Born  approximation. 


Modified  Born  Approximation  B 
This  is  nearly  the  same  as  Born  approximation  4  of 
Ref.  15.  We  no  longer  separate  f(0)  into  f  (0)  and  f.(0) 

3,  D 

factors,  but  instead  calculate  one  scattering  amplitude  fi 
for  r  <(  Rq  and  another  f^  for  r  )>  Rq,  and  then  combine  the 
two  coherently:  f  =  fi  +  fg.  In  f.,  q.^  again  replaces  qQ  (as 
opposed  to  the  procedure  of  Ref.  15)  and  the  nuclear  attenuation 
factor  also  appears.  Thus 
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However  in  fQ,  in  which  only  the  Coulomb  potential  appears,  the 
nuclear  attenuation  factor  is  dropped  and  the  outside  q  =  is 
used, 

f „  -  2?  /"  slnV  Ze2  a 

0  K2  R  -  -  dr 

“  o  q0r  r 

Choosing  a  nuclear  potential  V  =(-30  -  i  20)Mev  leads  to  the  dotted 
curve  in  Figure  8.  By  interpolating  between  the  curves  of  Figure  8, 
it  is  easily  seen  that  the  two  modified  Born  approximations  agree 
fairly  well,  and,  with  regard  to  angular  locations  of  maxima  and 
minima,  bofh  are  fairly  good  approximations  to  the  exact  phase 
shift  results. 


VII.  CONCLUSIONS 

Taking  into  account  the  limited  experimental  angular  and 
energy  resolution,  reasonably  good  agreement  is  obtained  for 
angles  up  to  80°  between  the  experimental  distribution  and  calcu¬ 
lated  curves  for  V  =  (-45  -  i20)  Mev  and  V  =  (-35  -  i20)  Mev.  The 


calculated  curves  for  V  =  (-45  -  i20)  Mev  give  slightly  better 
agreement  with  the  magnitude  of  in  the  forward  direction 

and  the  angular  location  of  the  first  diffraction  maximum. 


However,  the  calculated  results  are  expected  to  be  sensitive  to 

"I  /O  *1  O 

the  choice  of  nuclear  radius,  chosen  here  as  =  1.414  A  '  x  10  cm, 
and  the  shape  of  the  nuclear  density  distribution.  Introducing 
a  discontinuity  in  the  logarithmic  derivative  matching  at  the 
nuclear  surface  has  a  pronounced  effect  upon  the  shape  of  the 
calculated  curves,  but  n  values  less  than  about  0.75  do  not 
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appear  to  improve  the  agreement  "with  the  experimental  distri¬ 
bution,  at  least  "when  applied  to  the  square  well  model.  At 
backward  angles  where  incoherent  transitions  may  mask  the 
elastic  scattering,  the  predicted  maxima  and  minima  do  not 
appear  experimentally.  It  is  expected  that  a  more  reasonable 
nuclear  model  with  a  tapered  edge,  and  an  improvement  in  the 
experimental  angle  and  energy  resolution,  may  account  for  a 
large  part  of  the  discrepancy  that  still  exists. 
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FIGURE  CAPTIONS 

1.  Experimental  arrangement. 

2,  Experimental  angular  distribution  for  78  Mev  v~  and  v+ 
scattering  from  copper. 

3a.  Plot  of  2k^f  =  q  +  ip  for  optical  model  phase  shift  analysis 
of  v~  scattering  from  9  -  20°  to  80°  where  j f  j 2  =  ^  .  Curves 
are  for  =  -35  Mev,  V2  =  -20  Mev,  and  for  wave  function 
logarithmic  derivative  matching  factors  of  n  =  0.5,  0.75, 
and  1.00. 

3b.  Plot  of  2kgf  =  q  +  ip  for  9  =  40°  to  160°,  V,  =  -35  Mev, 

V2  =  -20  mev,  and  n  =  0.50,  0.75,  1.00. 

4.  Optical  model  solutions:  angular  distribution  for  77  Mev 
ir~  mesons  scattered  by  copper  for  =  -35  Mev,  V2  =  -20 
Mev.  a)  n  =  1.00;  b)  n  =  0.75;  c)  n  =  0.50. 

5a.  Plot  of  2kgf  =  q  +  ip  for  9  =  20°  to  80°,  =  -45  Mev, 

V2  =  -20  Mev,  and  n  =  0.50,  0.75,  1.00. 

5b,  Plot  of  2k^f  -  q  +  ip  for  9  =  40°  to  160°,  V-^  =  -45  Mev, 

V2  =  -20  Mev,  and  n  -  0.50,  0.75,  1.00. 

6.  Optical  Model  solutions:  angular  distribution  for  77  Mev 

ir~  mesons  scattered  by  copper  for  =  -45  Mev,  V2  =  -20  Mev. 
a)  n  =  1.00;  b)  n  =  0.75;  c)  n  =  0.50. 
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FIGURE  CAPTIONS  -  continued 

7.  Optical  model  solutions:  angular  distribution  for  78  Mev 
ir~  mesons  scattered  from  copper,  for: 

a)  V1  =  -45  Mev,  V2  =  -20  Mev,  n  =  1 

b)  VJL  =  -35  Mev,  V2  =  -20  Mev,  n  -  1 

c)  V-l  =  -25  Mev,  V2  =  -20  Mev,  n  =  1 

d)  Vx  =  -25  Mev,  V2  =  -15  Mev,  n  =  1 

8.  Modified  Born  approximation  results:  angular  distribution 
of  78  Mev  k~  mesons  scattered  by  copper,  for: 

a)  =  -45  Mev,  V2  =  -20  Mev,  Modified  Born  approximation  A 

b)  Vj  =  -35  Mev,  V2  -  -20  Mev,  Modified  Born  approximation  A 

c)  =  -25  Mev,  V2  =  -20  Mev,  Modified  Born  approximation  A 

d)  V1  =  -30  Mev,  V2  =  -20  Mev,  Modified  Born  approximation  B 
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TABLE  CAPTIONS 

I.  Experimental  angular  distribution  of  78  Mev  ir- 
scattered  from  copper. 

+ 

II.  Experimental  angular  distribution  of  78  Mev  v 
mesons  from  copper. 

III.  Calculated  optical  model  reaction  (total  minus 
cross  sections  for  78  Mev: 

=  kx  +  ik2  ,  V  =  V  +  iV2 

R0  =  5.658  x  10-13  cm 

k0RQ  =  4.80 

irR02  =  1006  mb. 


mesons 


scattered 


coherent) 
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TABLE  I 


glab 

20  +  4.8° 
25.25  +  4.9° 
30  +  4.6° 

34.5  +  3.9° 
40+3.9° 

44.5  +  4.2° 
47  +  4.2° 
50  +  4.3° 
55  +  4.3° 

57.5  +  5.8° 
60  +  4.4° 
70  +  5.8° 
75  +  5.4° 
77  t  6-£’° 
80  +  5.1° 
89  +  5.3° 

100  +  5.7° 
110  +  6.0° 
120  +5.9° 
130  +  6.4° 
145  +  4.4° 


jgr  (mb) 
1459  +  89 

747  +  78 

385  +  31 

177  +  27 

59.4  +  5.5 

57.5  +  13.2 

41.4  +  5.0 

36.1  +  3.6 

33.1  +3.7 
31.3  +  4.5 
27.7  +  1.6 

12.12  +  0.89 
11.52  +  1.18 
11.74  +  2.37 
12.01  +0.81 
11.58  +  0.84 
10.39  +  0.84 
9.13  +  0.74 
10.19  +  1.46 
0 • 95  +  1 . 14 
11.04  +  1.41 
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TABLE  II 


e 


lab 

m  (rab> 

15  +  2.9° 

1882  +  342 

20  +  5.0° 

748  +  96 

30  +  4.2° 

267  +  44 

40  +  4.2° 

63.3  +  13.4 

50  +  4.3° 

* — 

30.2  +  8.6 

60  +  5.9° 

16.2  +  6.1 

HO  +  6.5° 

8.29  +  0.84 

TABLE  III 


(567) 


Fig.  1.  Experimental  Arrangement 


